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ABSTRACT 

We provide a model for how Kennicutt-Schmidt (KS) laws, which describe the correlation between 
star formation rate and gas surface or volume density, depend on the molecular line chosen to trace 
the gas. We show that, for lines that can be excited at low temperatures, the KS law depends on how 
the line critical density compares to the median density in a galaxy's star-forming molecular clouds. 
High critical density lines trace regions with similar physical properties across galaxy types, and this 
produces a linear correlation between line luminosity and star formation rate. Low critical density 
lines probe regions whose properties vary across galaxies, leading to a star formation rate that varies 
superlinearly with line luminosity. We show that a simple model in which molecular clouds are treated 
as isothermal and homogenous can quantitatively reproduce the observed correlations between galactic 
luminosities in far infrared and in the C0(1 0) and HCN(1 0) lines, and naturally explains why 
these correlations have different slopes. We predict that IR-line luminosity correlations should change 
slope for galaxies in which the median density is close to the line critical density. This prediction 
may be tested by observations of lines such as HC0+(1 0) with intermediate critical densities, or 
by HCN(1 0) observations of intensely star-forming high redshift galaxies with very high densities. 
Recent observations by Gao et al. hint at just such a change in slope. We argue that deviations from 
linearity in the HCN(1 0)— IR correlation at high luminosity are consistent with the assumption of 
a constant star formation efficiency. 

Subject headings: ISM: clouds — ISM: molecules — stars: formation — galaxies: ISM — radio lines: 
ISM 
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strong correlation between the IR luminosity of galax- 
ies and emission in the HCN(1 0) line, which mea- 
sures the mass at densities significantly greater than 
that probed by C0(1 — > 0). However, they find that 
their correlation, which covers nearly three decades in 
total g alactic sta r formation rate, is linear: Lfir oc 
Lhcn- IWu et al] ([2005) show that this correlation ex- 
tends down to individual star-forming clumps of gas in 
the Milky Way, provided that their infrared luminosities 
are ^ 10'* ''^ Lq. Interestingly, however. iGao et al] ()2007f ) 
find a deviation from linearity in the IR-HCN correlation 
for a sample of intensely star-forming high redshift galax- 
ies. These sources show small but significant excesses of 
infrared emission for their observed HCN emission. 

The difference in power law indices between the Lfir — 
Leo a-nd LpiR — Lhcn correlations is statistically signif- 
icant, and, on its face, puzzling. An index near N = 1.5 
seems natural if one supposes that a roughly constant 
fraction of the gas present in molecular clouds will be 
converted into stars each free-fall time. In this case one 
expects cx pj-^ (jMadord fl977l : lElmegreenlll994f ). If 
gas scale heights do not vary strongly from galaxy to 
galaxy, this implies S* oc S^'^ as well, which is consis- 
tent with the observed Kennicutt law. More generally, 
since the dynamical timescale in a marginally Toomre- 
stable (Q ~ 1; see iMartin fc Kennicuttj l2001h galactic 
disk is of order fi"^ cx (Gpg)"^/^, where fl is the angular 
frequency of the disk, an index close to = 1.5 is ex- 
pected if star formation is regulated by any phenomenon 
that converts a fixed fraction of the gas into stars on this 

Electronic address: |krumholz@astro.princeton.edu, thompQastro.princettin gditCale (jElmegrecn 2002)^ 

2 Hubble Fellow Qn the Other hand, Wu et~all ()2005D suggest a simple 

Lyman Spitzer Jr. Fellow interpretation of the linear IR-HCN correlation. They ar- 



1. INTRODUCTION 

ISchmidtj ()1959l Il963f) first proposed that the rate at 
which a gas forms stars might follow a simple power 
law correlation of the form oc , where is the 
star formation rate per unit volume, pg is the gas den- 
sity, and N is generally taken to be in the range 1 — 2. 
In the decades since, observations have revealed two 
strong correlations that appear to be evidence for this 
hypothesis. First, galaxy surveys reveal that the in- 
frared luminosity of a galaxy, which traces the star for- 
mation rate, varies with its luminosity in the C0(1 —^ 0) 
line, which trace s the total mass of mo lec ular gas, as 
Lftr oc Lhj-^ '^ (IGao k SolomonI l2004alfbl: iGreve et all 
120051: iRiechers et al.ll2006a[ ). iKennicuttI (|1998allbl ) iden- 
tified the closely-related correlation between gas surface 
density Eg and star formation rate surface density E,, 
E* oc Yig-'^^'^-^^ , a relation that has come to be known 
as the Kennicutt Law. Since over the bulk of the dy- 
namic range of Kennicutt 's data galaxies are predomi- 
nantly molecular, this is effectively a correlation between 
molecular gas, as traced by C0(1 0) line emission, 
and star formation. Spatially resolved observations of 
galaxies confirm that, at least for molecule-rich galax- 
ies where resolved C0(1 — > 0) observations are possible, 
star formation is more closely coupl ed with ga s traced by 
C0(1^ 0) than with atomic gas (IWongfc B litz 2002; 
Hever et al.l l2004 iKomugi et all l2005t lKe°icutt et al.i 
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Second, iGao fc SolomonI ()2004al lbll find that there is a 
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gue that the individual HCN-emitting molecular clumps 
that they identify in the Milky Way represent a funda- 
mental unit of star formation. The linear correlation 
between star formation rate and HCN luminosity across 
galaxies arises because a measurement of the HCN lu- 
minosity for a galaxy simply counts the number of such 
structures present within it, each of which forms stars at 
some roughly fixed rate regardless of its galactic environ- 
ment. However, in this interpretation it is unclear why 
the structures traced by HCN(1 0) emission should 
form stars at the same rate in any galaxy. After all, one 
could equally well argue that molecular clouds traced by 
C0(1 0) are fundamental units of star formation, but 
the non-linear IR-CO correlation clearly shows that these 
objects do not form stars at a fixed r ate per uni t mass. 
Moreover, the evidence presented by iGao et all ()2007f ) 
that the linear IR-HCN correlation varies in extremely 
luminous high redshift galaxies suggests that the rela- 
tionship between HCN emission and star-formation may 
be somewhat more complex. 

In this paper we attempt to explain the origin of the 
difference in slope between the CO and HCN correla- 
tions with star formation rate, and more generally to 
give a theoretical framework for understanding how cor- 
relations between star formation rate and line luminosity, 
which we generically refer to as Kennicutt- Schmidt (KS) 
laws, depend on the tracer used to define them. Our cen- 
tral argument is conceptually quite simple, and in some 
sense represents a combination of the intuitive arguments 
for CO and HCN given above. 

Consider an observation of a galaxy in a molecular 
tracer with critical density ricrit, which essentially mea- 
sures the mass of gas at densities of ricrit or more, i.e. the 
gas that is dense enough for that particular transition to 
be excited. In galaxies where the median density of the 
molecular gas is significantly larger than ncrit, this means 
that the observation will detect the majority of the gas, 
and the bulk of the emission will come from gas whose 
density is near the median density. Since the gas den- 
sity will vary from galaxy to galaxy, the star formation 
rate per unit gas mass will vary as roughly p^'^ , with one 
factor of pg coming from the amount of gas available for 
star formation, and an additional factor of p^'^ coming 
from the dependence of the free-fall or dynamical time 
on the density. 

On the other hand, in galaxies where the median gas 
density is small compared to the critical density for the 
chosen transition, observations will pick out only high 
density peaks. Since the density in these peaks is set by 
ricrit, and not by the conditions in the galaxy, these peaks 
are at essentially the same density in any galaxy where 
they are observed, and the corresponding free-fall times 
in these regions are constant as well. As a result, the star 
formation rate per unit mass of gas traced by that line 
is approximately the same in every galaxy, because the 
corresponding free-fall time is the same in every galaxy. 

In the rest of this paper, we give a quantitative version 
of this intuitive argument, and then discuss its conse- 
quences. In § [2] we develop a simple formalism to com- 
pute the star formation rate and the molecular line lu- 
minosity of galaxies, and in §[3] we use this formalism to 
predict the correlation between star formation rate and 
luminosity. We show that our predictions provide a very 



good fit for a variety of observations, and make predic- 
tions for future observations. We discuss the implications 
of our work and its limitations in § |31 and summarize our 
conclusions is § [5] 

2. STAR FORMATION RATES AND LINE LUMINOSITIES 

2.1. Cloud Properties 

Consider a galaxy in which the star-forming molecu- 
lar clouds have a volume-averaged mean molecular hy- 
drogen number density n = Jj /^h2j where p is the 
volume-averaged mass density of the molecular clouds in 
the galaxy and fj,H2 = 3.9 x 10"^"' g is the mean mass per 
hydrogen molecule for a gas of standard cosmic composi- 
tion. Observations indicate that n varies by two to three 
decades over the galaxies for which the Kennicutt and 
Gao & Solomon correlations are measured, from n » 50 
cm"^ in normal spirals like the Milky Way (jMcKed 
Il999| ) up to n « 10* cm~'^ in the strongest starburst 
syste ms in the local universe (e.g. [Downcs & Solomo^ 
|1998[ ). There is strong evidence that densities in molec- 
ular clouds follow a lognorm al probability distribution 
fun ction (PDF; see reviews by iMac Low fc Klessen|[200l 
and lElmegreen &: ScaL 
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where a; = n/n is the molecular hydrogen number den- 
sity n relative to the average density, a is the width of 
the lognormal, and Ina; = For this distribution 

the median density is rimed = n &cp{(j^ / 2) . Numerical 
experiments show that for supersonic isothermal turbu- 
lence cr^ ~ In (1 + 3A4V4), where M is the ID Mach 
number of the turbulence (Nordlund fc Padoan' 1999'; 
lOstriker et"an 119991 : [Padoan fc Nordlund 2002). Mach 
numbers i n star-formin g molecular clouds range from 
- 30 (lMcKeel ll999D in normal spirals to M - 100 
in strong starbursts (jPownes fc Solomon! Il998h . imply- 
ing that median densities in molecular clouds range from 
^ lO'' cni^"^ in normal spirals to ~ 10^ cm"'^ in star- 
bursts. Star forming clouds within a galaxy are approx- 
imately isothermal, except very near strong sources of 
stellar radiation, so we assume a fixed temperature T 
for the clouds. Observati onally, T ran ges from roughly 
10 K in normal spirals (jMcKed 119991) up to as much 
about 50 K in stroiig starb ursts (jDownes fc SolomonI 
[19981: iGaofc Solomonll2004bl) . 

2.2. Star Formation Rates 

First le t us ask how quick l y star s form in such a 
medium. iKrumholz fc McKed (|2005[ ) give a model for 
star formation regulated by supersonic turbulence in 
which a population of molecular clouds of total mass Md 
form stars at a rate M* = SFRg Md/tg (n), where iff (n) 
is the free-fall time evaluated at the mean density and 
SFRff is a number of order 10"^ that depends weakly on 
M.. We therefore estimate the star formation rate per 
unit volume as a function of the mean density given by 



SFRff 1 



(2) 



We adopt the Krumholz fc McKee result SFRff sa 
0.014(Al/100)""-^2 fQj. clouds with a fiducial virial ra- 
tio of Q!vir = 1-3. 
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Alternately, iKrumholz fc TanI ()2007f ) point out that 
observed correlations between the star formation rate 
and the luminosity in different density tracers imply that 
over a 3 — 4 decade range in density n, 



w 10 



_2 Afci(> n) 
ts{n) 



(3) 



where Mci(> n) is the mass of gas with a density of n 
or higher, and Md = Mc\{> 0). For a given choice of 
n this provides an alternative estimate of the star for- 
mation rate which is purely empirical, and independent 
of any particular theoretical model. However, the differ- 
ence between the star formation rates predicted by ([2]) 
and ([3|) is small. For gas with a lognormal PDF, 



Mei(> n) 
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1 -I- erf 
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and using this to evaluate equation ([3]) indicates that, for 
Mach numbers in the observed range, the two prescrip- 
tions ([2]) and ([3]) give about the same star formation rate 
over a very broad range in x. For example, at = 30 
the two estimates agree to within a factor of 3 for den- 
sities in the range 0.2 < a; < 4 x 10''. Given the scatter 
inherent in observational estimates of the star formation 
rate, a factor of 3 difference is not particularly signifi- 
cant, so it matters little which prescription we adopt. In 
practice, we will use equation ([2]). 

2.3. Line Luminosities 

Now we must compute the luminosity of molecular line 
emission from the galaxy. Even for a cloud that is not 
in local thermodynamic equilibrium (LTE), for optically 
thin emission this calculation is straightforward. How- 
ever, the molecular lines used most often in galaxy sur- 
veys are generally optically thick. To handle the effect of 
finite optical depth on molecule level populations and line 
luminosities, we adopt an escape probability approxima- 
tion and treat clouds as homogeneous spheres. This is 
not fully consistent with our assumption that clouds have 
lognormal density PDFs, since the escape probability for- 
malism assumes a uniform level population throughout 
the cloud, and the essence of our argument in this paper 
turns on how the level population varies with density. 
However, this approach gives us an approximate way of 
incorporating the optical thickness of star-forming clouds 
into our model, the only alternative to whi ch for turbu- 
lent m edia is full numerical simulation (e.g. lJuvela et al.l 
l2001h . We therefore proceed by treating clouds as ho- 
mogeneous in order to determine their escape probabili- 
ties, and we then relax the assumption of homogeneity, 
while keeping the escape probabilities fixed, in order to 
determine level populations and cloud luminosities as a 
function of density. 

Consider a cloud of radius R in statistical equilibrium 
but not necessarily in LTE. In the escape probability 
approximation, the fraction fi of molecules of species S 
in state i is given implicitly by the linear system 
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h (5) 
(6) 



where qij is the collision rate for transitions from state 
i to state j, Aij is the Einstein spontaneous emission 
coefficient for this transition, /3y is the cloud-averaged 
escape probability for photons emitted in this transition, 
the sums are over all quantum states, and we understand 
that Aij — for i < j and qij — for i — j. 

Equations (O and (O allow us to compute the level 
populations fi for given values of (3ij. To completely 
specify the system, we must add an additional consis- 
tency condition relating the values of f3ij to the level 
populations. For a homogeneous spherical cloud, the es- 
cape probability for a given line is related to the optical 
depth from the center to the edge of the cloud by 
(B. Draine, 2007, private communication) 
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1 + 0.5t,, 



(7) 



where Tij is computed at the central frequency of the 
line. In turn, the optical depth is related to the level 
populations by 



9j 
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5j 4(2^)3/2 A^c, V I39^ 



(8) 



where Xij is the wavelength of transition i j, gi and 
gj are the statistical weights of states i and j, Cs is 
the isothermal sound speed of the gas, and X{S) is the 
abundance of molecules of species S. Note that this 
equation implicitly assumes that the cloud has a uni- 
form MaxwcUian velocity distribution with ID disper- 
sion Mcs, consistent with our treatment of the clouds 
as homogeneous spheres. One additional complication is 
that we do not directly know cloud radii for most exter- 
nal galaxies, where observations cannot resolve individ- 
ual molecular clouds. However, we often can diagnose 
the optical depths of transitions by comparing line ratios 
of molecular isotopomers of different abundances. We 
therefore take tiq , the optical depth of the transition be- 
tween the first excited state and the ground state, as 
known. For a given level population this fixes the value 
of R. 

We solve equations ©-([8]) using Newton-Raphson it- 
eration. In this procedure, we guess an initial set of es- 
cape probabilities (3ij , and solve the linear system ^ and 
([6]) to find the corresponding initial level populations fi . 
We then compute the optical depths from equation 
([8|). The guessed escape probabilities Pij and the corre- 
sponding optical depths generally will not satisfy the 
consistency condition ([7]) , so we then iterate over /3y val- 
ues using a Newton-Raphson approach, seeking Pij for 
which the level populations give optical depths Tij such 
that all elements of the matrix Pij — 1/(1 + O.St^) are 
equal to zero within some specified tolerance. We use 
the LTE level populations and escape probabilities for 
our initial guess, so that the iteration converges rapidly 
when the system is close to LTE. 

Once we have determined the escape probabilities (3ij , 
we compute the luminosity by holding the Pij values fixed 
but allowing the level populations to vary with density, 
then integrating over the PDF. Thus, the total luminos- 
ity per unit volume in a particular line is 



L,j = X{S)l3,jAjh 



fin 



dp 
d\na 



dlnx, 



(9) 
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TABLE 1 

Model Parameters 



Parameter 


Normal galaxy 


Intermediate 


Starburst 


Reference 


T 


10 


20 


50 


1-4 


M 


30 


50 


80 


1-4 


X(CO) 


2 X 10-4 


4 X 10-1 


8 X lO-'' 


5 


X(HCO+) 


2 X 10-9 


4 X 10-9 


8 X 10-9 


6, 7 


X(HCN) 


1 X 10-* 


2 X 10-** 


4 X 10-* 


6-8 


'^CO{1-.0) 


10 


20 


40 


9 


■^HCO+(l-.0) 


0.5 


1.0 


2.0 


6, 7 


■'"HCN{1^0) 


0.5 


1.0 


2.0 


6, 7 


OPR 


0.25 


0.25 


0.25 


10 



Note. — OPR = H2 ortho- to para-ratio. References: 
1 - Solomon et al. (1983), 2 - [Uao fc Solomon (2004b), 3 - 
Downes fc Solomon (199 ^4 - | Wu et al.l ( 2005) , 5 - Black (2000), 6 - 
Nguyen et al. (1992). 7 - Wild et al.l (119921). 8 -ILaliuis fc wm DishoecM 



(120001) . 9 - Combes (1991'). 10 - INeufeld et"an (I2OO6I) 



where Vij is the hne frequency, fi is an impHcit function 
of n given by the solution to equations ([5|) and ([6]), and 
we assume that the abundance X(S) is independent of 
n. The line luminosity per unit mass is /(/i/z^ri). 

An IDL code that implements this calcu- 
lation is available for public downlo ad from 
http:/ /www. astro. princeton.edu/~krumholz/ astron- 
omy, html. 

3. CORRELATIONS AND KENNICUTT-SCHMIDT LAWS 
3.1. Lines and Parameters 

Using the formalism of § [21 we can now predict the 
correlation between the star formation rate and the 
luminosity of a galaxy in molecular lines. We make 
these predictions for three representative molecular lines: 
C0(1 ^ 0), HC0+(1 ^ 0), and HCN(1 ^ 0). For 
the first and last of these transitions, there are exten- 
sive observational surveys. We select HC0+(1 0) 
in addition to these two because there is some obser- 
vational data for it, and because its critical density of 
= /3hco+4.6 X 10 cm makes it intermediate 
between C0(1 — > 0), with ricrit = Pco^QO cm~'^, and 
HCN(1 0), with Merit = /3hcn2.8 x 10^ cm'^.s Here 
Ps is the escape probability for the 1^0 transition of 
species S. These critical densities are for T = 20 K. All 
molecular data are t aken from the Leiden Atomic and 
Molecular Database'* (jSchoier et al.| [2005). 

We make our calculations for three sets of fiducial pa- 
rameters which we summarize in Table [T] The three sets 
correspond roughly to typical conditions in normal disk 
galaxies like the Milky Way, to starburst galaxies like 
Arp 220, and to a case intermediate between the two. We 
have selected parameters for each case to roughly model 
the systematic variation of ISM parameters as one moves 
from normal disk galaxies to starbursts. Thus, we vary 
the ISM temperature from 10 — 50 K and the molecular 
cloud Mach number from 30 — 80 as we move from Milky 
Way-like molecular clouds to tem peratures and Mach 
numbers typical of starbursts (e.g. iDownes fc Solomonl 

^ Note that our critical dens ity for HCNfl — > 0) is somewhat 
larger than the value quoted bv lGao &: Solomonl (I2004a irbl) . proba- 
bly because their calculation is based on somewhat different as- 
sumptions about how to extrapolate from calculated rate coef- 
ficients for HCN c ollisions with He to collisions with H2. See 
ISchoier et aP ((20051 ) for details. 

* http:/ /www. strw.leidenuniv.nl/~moldata/ 



|1998[) . Similarly, starbursts, which preferentially occur 
at galactic centers, have systematically l arger metallici- 
ties than galaxies like the Milky Way (e. g. IZaritskv et al.l 
ll994HYao et"aIll2003l : lNetzer et al.ll2005[ ). To explore this 
effect, we use abundances and 1 — > optical depths are 
twice and four times as large for our intermediate and 
starburst models, respectively, as for our normal galaxy 
model. 

3.2. Kennicutt- Schmidt Laws 

We first plot, in Figure [H the quantities L^^[dL{< 
n)/d\'Q-n] [solid lines) and M~^[dM{< ri)/d\a.n] [dot- 
ted lines) as a function of density n for galaxies with 
mean densities n = lO'^, lO'^, and 10'' cm~'^, for the 
tracers C0(1 ^ 0), HC0+(1 ^ 0), and HCN(1 ^ 
0), and for the Mach number and temperature corre- 
sponding to our intermediate case in Table [TJ Here 
L[< n) and M[< n) are the luminosity and mass per 
unit volume contributed by gas of density n or less, 
i.e. L[< n) = X[S)PijAijhi'ij J^_^^ fin[dp/dliin)d\iin, 
M[< n) = fiH^n[dp/dlnn)d\iin, L = L[< 00), and 

M ^ M[< 00). Physically, L"i[dL(< n)/d\nn] and 
AI^^[dM[< n)/dlnn] represent the fractional contribu- 
tion to the total line luminosity and the total mass that 
comes from each unit interval in the logarithm of den- 
sity. The plot shows what density range provides the 
dominant contribution to the line luminosity in differ- 
ent lines and for galaxies of differing mean densities, and 
how the gas contributing light compares to the gas con- 
tributing mass. Because the mass distribution is entirely 
specified by n and Ai, the dotted lines are the same in 
each of the three panels. Additionally, because of our 
choice = 50 (Table [T|), the median density (the den- 
sity corresponding to the peak in M~^[dAI[< n)/dlnn]) 
is ".mod ~ 4371. In each panel, the critical density for 
each molecule is identified by a vertical dashed line. 

The top panel clearly shows that for the CO line, the 
light and the mass track one another very closely, even 
at the lowest densities. Thus, because n,ncd > JT-crit, the 
solid lines move in lock-step with the dashed lines as n 
increases. In contrast, for HCN most of the luminosity 
comes from densities near the critical density regardless 
of the mass distribution. For the lowest n this means that 
the line luminosity is entirely dominated by the high den- 
sity tail of the mass distribution. As the median density 
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Fig. 1. — Fractional contribution to the total luminosity 
L~^[dL{< n)/dlnn] (solid lines) and mass M~^[dM{< n)/dlnn] 
(dotted lines) versus density n for the lines CO(l 0) (top 
panel), HCO+(l 0) (middle panel), and HCN{1 0) 

(bottom panel). The three curves show the cases n = 10^ 
cm~^, 10^ cm~^, and 10"* cm"'', from leftmost to rightmost. 
We also show the critical density of each molecule, corrected 
for radiative trapping (dashed vertical lines). These calcula- 
tions use the parameters for the intermediate case listed in Table [l] 

rimed varies by a factor of 100 (from 4.3 x 10^ — 4.3 x 10^ 
cm~^), the peak of L~^[dL(^< n)/dlnn] moves by just 
a factor of a few in n. The HCO+ line is intermedi- 
ate between CO and HCN. For n = 10^ cm~^ and 10^ 
cm~^, rimed ~ '^crit, and as with HCN most of the emis- 
sion comes from near the critical density. For n = 10* 
cm~^, rimed > ricrit, and the light starts to follow the 
mass, in a pattern similar to that for CO. Although Fig- 
ure[T]shows only the intermediate case, the normal galaxy 
and starburst cases give qualitatively identical results. 
This confirms the intuitive argument given in § [TJ high 
critical density transitions trace regions of similar den- 
sity in every galaxy, while low critical density transitions 
trace regions whose density is close to the median density. 

Now consider how the luminosity in a given line corre- 
lates with the star formation rate in galaxies of varying 
mean densities. For a given n, we can compute the vol- 
ume density of star formation from equation Q and the 
line luminosity density from To facilitate compari- 
son with observations, rather than consi dering the total 
line l uminosity, we use the quantity L' (jSolomon et al.l 
|1997[) . which is related to the luminosity L by 



V 



(10) 



converted to the units K km s^^ pc^. 

Similarly, we can estimate the far infrared luminosity 
from the star formation rate. There is a tight correla- 
tion between far-IR emission and star formation, par- 
ticularly for dense, dusty galaxies li ke those th a t make 
up most of the dynamic range of the lKennicuttI (|1998af) 



sample. To the extent that most or all of the light from 
young stars is re-processed by dust before escaping the 
galaxy, the bolometric luminosity integrated over the 
wavelength range 8 — 1000 /im, which we define as Lfir, 
simply provides a calorimetric measurement of the total 
energy output by young stars, and is therefore an excel- 
lent tracer of recent star formation (ISanders fc Mirab e ll 

1 QQfil- iTJnwan-RoKi'riQr.n ot al ' 'l QQTl 'Tf onm'mittl M QQSallVj- 



19961: iRowan-R obinson et al.' '1997": 'Kennicutt' '1998a 
Hirashita et al.l 2003; BcD 2003; Iglesias-Paramo et al.l 



2006f ). We therefore estimate the FIR luminosity from 



the star formation rate via 



^FIR 



(11) 



wher e e is an IMF-dep endent constant. For consistency 
with lKennicuti (|1998al |bh. we take e = 3.8 x 10""*. To be 



precise and to facilitate comparison with observations, we 
adopt the Sanders & Mirabel (1996) definition of L-pm, 
as a weighted sum of the luminosity in the 60 and 100 
/im IRAS bands. This definition of the infrared luminos- 
ity generically underestimates the total inf rared luminos- 
ity [8 - 1000] /iin by a factor of 1.5 - 2 (jCalzetti et al.l 
120001 : iDale et all 1200 ll : 03cll 2003). However, we use the 
e value appropriate for Lfir rather than for the total 
IR luminosity because some of the observations to which 
we wish to compare our model (see § 13. 3p provide only 
L-piK- Note that this choice for the connection between 
the star formation rate and the infrared luminosity is not 
fully consistent with our choice of the gas temperature 
for the three sets of parameters — normal, intermediate, 
and starburst — listed in Table [1] an issue we discuss in 
more detail in § 14.31 

We plot the ratio of star formation rate to line lumi- 
nosity, and infrared luminosity to line luminosity, as a 
function of n in Figure [2l First consider the top panel, 
which shows all three lines computed for the interme- 
diate case. This again confirms our intuitive argument. 
Since the luminosity per unit volume in the CO line is 
roughly proportional to the mass density, and the star 
formation rate / IR luminosity is proportional to mass 
density to the 1.5 power, the ratios M^jL' and L-pYs,/L' 
vary roughly as rfi^^ . A powerlaw fit to the data over the 
range shown in Figure [2] gives an index of 0.57. In con- 
trast, the ratio of star formation density to HCN luminos- 
ity density is nearly constant for galaxies with n < 10"^ 
cm^"^, and varies quite weakly with n up to densities 
of 10'* cm"'^, values found in the densest starbursts. A 
powerlaw fit from 10 cm^'^ to 10* cm~^ gives an index of 
0.17; from 10 cm~^ and 10'^ cm~'^, the best fit powerlaw 
index is 0.08. As in Figure [TJ the slope of the M*/L' 
curve for HCO^ represents an intermediate case, with 
a roughly constant ratio of A/*/L' and Lpi^/L' at low 
n, rising to a slope comparable to that for CO at high 
values of n. 

Now consider the bottom three panels in Figure [21 
Each panel shows the ratio of star formation rate and 
infrared luminosity to line luminosity for a single line, 
computed for each of the three galaxy models. The most 
important point to take from these plots is that the choice 
of galaxy model has little effect in most cases. The largest 
differences are for HCN, where at ri = 10 cm~'^ the IR to 
line ratio predicted for the intermediate case differs from 
the normal galaxy case by a factor of 6.1, and from the 
starburst case by a factor of 4.1. This variation comes 
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Fig. 2. — Ratio of star formation rate or infrared luminosity to 
line luminosity, as a function of mean density n. In the top panel 
we show the lines CO(l 0) (solid line), HCO+(l 0) [dot- 
dashed line), and HCN(1 — > 0) (dashed line) for the intermediate 
case in Table[T] In the next three panels we show the CO(l — > 0), 
HCO^(l — > 0), and HCN(1 — > 0) lines for the normal galaxy case 
(dot-dashed line), intermediate case (solid line), and starburst case 
(dashed line). 

primarily from changes in the Mach number and the op- 
tical depth between models. The higher Mach number of 
the starburst model significantly increases the amount of 
mass in the high ovcrdensity tail of the probability dis- 
tribution, while the higher optical depth lowers the ef- 
fective critical density. Both of these effect increase the 
amount of mass dense enough to emit in HCN(1 0) 
and reduce M^:/L' . At higher mean densities these effects 
become less important and the models converge, so that 
by n = 10"* cm~'^ the range in M*/L' from the normal 
to the starburst case is only a factor of 3.5. 

Most importantly, our central conclusion that 
M,/i{j(-;j^ is roughly constant across galaxies, while 
M,/L^Q rises as roughly [L'^q]'^-^, still holds when we 
consider how conditions vary across galaxies. Galaxies 
with low mean densities n are generally closest to the 
normal galaxy case, while those with high mean densi- 
ties should be closest to the starburst case, and this sys- 
tematic variation in galaxy properties with n still leaves 
M^/L' relatively flat for HCN, and varying with a slope 
close to 0.5 for CO. From the normal galaxy case at 
n = 10 cm"^ to the starburst case at n — 10"* cm~^, 
the value of M„/L' varies by more than a factor of 50 
for the C0(1 0), but by less than a factor of 3 for the 
HCN(1 -> 0). 



3.3. Comparison with Observations 

The calculations illustrated in Figure [2] demonstrate 
the basic argument that one expects a roughly constant 
star formation rate per unit line luminosity for high den- 
sity tracers (e.g., HCN), and a star formation rate per 
unit luminosity that rises like luminosity to the ~ 0.5 
for low density tracers (e.g., CO). However, in large sur- 
veys one cannot always determine the mean density in 
a galaxy, which would be required to construct an ob- 
servational analog to Figure [H Instead, we can use our 
calculated dependence of star formation rate and line 
luminosity on density to compare to observations as fol- 
lows. Equation ([9]) gives the total molecular line lumi- 
nosity per unit volume and equation ([2|) gives the star 
formation rate, which we convert to an IR luminosity 
via equation (jlip . For fixed assumed volume of molecu- 
lar star-forming gas (Knoi) we can then predict the ex- 
pected correlations between L' in a given molecular line 
and LpiR. The three panels of Figure [3] show our results 
for XpiR as a function of L'qq, ^hcn' ^'^'^ -^hco+ ^"-"^ 
intermediate model (see Table [IJ and for several values 
of Vinoi- Figure [4] shows how are results vary as a func- 
tion of the assumed T and M.. There, for fixed Knoi, 
we compare our predictions for the intermediate model 
with the normal and starburst models. In both figures we 
compare our models to dat a culled from the literature. 

From th e wo rk of iGao k, SolomonI f;2 004a,li') 
iGreve et all (|2005l t heir Fig. 7 ) iRiechers et"al.l (|2006bl . 
their Fig. 5), and iGao et all (|2007| ). as well as the 
theoretical arguments in the preceding sections, we 
expect a strong, but not linear, correlation between 
the CO luminosity and the star formation rate — 
as measured by Lfir — with the approximate form 
LpiR oc L'^q . The left panel of Figure [3] shows the 
CO data, the approximate correlation expected {solid 
line segment; offset from the data for clarity) and the 
theoretical prediction (solid lines) for a total volume of 
molecular gas of Vmoi = 10'', 10^, and 10^ pc'^. Because 
at fixed Lfir, galaxies exhibit a dispersion in Vmoi we 
expect there to be intrinsic scatter in this correlation, 
roughly bracketed by the range of Knoi plotted. 

The middle and right panels of Figure [3] show the same 
prediction for ^^00+ ^hcn- these cases, because 
the molecular line luminosity is nearly linearly propor- 
tional to LpiR, the dependence on V^noi is much weaker 
than for L'qq- However, systematic changes or differences 
in the fiducial parameters for the calculation (see Table 
[1]) introduce uncertainty and scatter into the correlation. 
Figure m assesses this dependence. It compares the pre- 
dictions of our model for normal {dot- dashed lines), inter- 
mediate {solid lines), and starburst {dashed lines) galax- 
ies, as defined in Table [H for fixed V^oi = lO'^ pc'^. Our 
simple model reproduces the data rather well, and it pre- 
dicts that generically there may be more intrinsic scatter 
in the -/jqq — Lfir correlation than in either ^hcn ~ ^fir 

Note that in both the middle and right panels of Fig- 
ures [3] and m one expects a turn upward in the corre- 
lation at high LpiR, a- deviation from linearity. This 
follows from the fact that in our model, at fixed Knoi, 
systems with higher Lfir have higher average densities. 
At sufficiently high Lfir we thus expect L^qj^— LpiR and 
L^pQ+— Lfir to steepen, in analogy with the L'qq — Lfir 




Fig. 3. — i/FlR {Lq) versus LJ-.q(1 — > 0) {left panel), ^^^^^^(1 — > 0) {middle panel), and ^^^^(1 — > 0) (K km s ^ pc^; right panel). 

The lines in each panel derive from the model presented in this paper with a constant total volume of molecular material of 10^, 10*, and 
10^ pc^ (lowest to highest). The thick solid line segmen t shows power-law slopes to guide the eye. Data in the left and right panels are 
from lGao fc SolomonI ll2004allB ( circles) and I Gao etall ||2007, ) ( open squares for detectio ns, arrows for upper limits). The midd l e pane l 
combines data from lNg uyen et al. (1992) {small circles with lines), Gracia-C arpio et al.l (2006) {big circles), and Ricch eFs et al.l Il2006bl ) 
{open square; using the IGao et al.. ( 2007) FIR luminosity and magnification factor). For all data, Lfir is define d based on a weighte d sum 
of the galaxy luminosity in the 60 and 100 /xm IRAS bands, as described by Sanders & Mirabel (1996). For the lNguven et all ||1992!) data, 
the uncertainties in Ljj(-.q+ indicated by the lines arise because i Nguyen ct al. provide both HCN(1 — > 0) a nd HC0"*"(1 — > 0) intensities, 
but the values for L'^^j^ derived from t heir work gene rally fall a factor of 2 — 3 below the i^j^jj from Gao fc Solomon ( 2004a, h) for the 
same systems. This is probably because INguven et al.l use a single beam pointing rather than integrat ing fully over e xtended sources, and 
therefore miss some of the flux. We therefore show two values of L^^^^ , connected by a line, for each INguven et al.l data point: a smaller 
value calculated directly from the data listed in their Table 2, and a larger value obtained by multiplying the L'^q^ value of lGao &: SolomonI 
for that galaxy by the ratio /hco+/^hcn measured by INguven et aLl lf this ratio is constant over the source, this estimate should correctly 
account for the flux outside the beam in the lNguven et all HCO+ observation. 



correlation. The data points with very high Lfir, in Fig- 
ures [3] and [H which might be used to test this prediction 
of our model, are gravitationally lensed, at high redshift, 
and contaminated by bright AGN. It is therefore un- 
clear if the deviation from linearity implied particularly 
by the upper limits in L^^j^ in the right panels of Fig- 
ures [3] and 2] is a result of enhance d Lfir, caused by the 
AGN emission (ICarim et al.li2005D . or is instead a result 
of less molecular line emission per unit st a r form ation, 
as our model implies (Fig. IGao et all (|2007[ ) note, 
however, that in the three systems for which the contri- 
bution from the AGN has been estimated (F10214-I-4724, 
D. D ownes fc P. Solom on 2007, in preparati on; Clover- 
leaf, iWei fi et al.l [200l APM 08279-^5255, iWeifi et all 
[2005, 200f the corrections are only significant for APM 
08279-1-5255. This suggests that the data are so far con- 
sistent with our interpretation, but clearly much more 
data at high Lfir — or, more precisely for our purposes, 
at high density — is required to test our predictions. We 
discuss the issue of AGN contamination further in § 14.31 
As a final note, the data so far do sup port the utility 
of HC O+ as a useful tracer of dense gas. iPapadopoulosI 
()2007f ) has argued against the utility of HCO+ as a faith- 
ful tracer of mass in starbursts on the basis that, since it 
is an ion, its abundance is strongly dependent on the free- 
electron abundance and might therefore vary strongly 
between galaxies with different ionizing radiation back- 
grounds. We cannot rule out this possibility given the 
limited data set available, but we see no strong evi- 
dence in favor of it from the data shown in Figures [3] 
and m As we have argued, HC0"'"(1 0) is particu- 
larly useful because its critical density is between that 
of C0(1 ^ 0) and HCN(1 0) and, thus, as Fig- 
ure [3] and [4] show, the correlation between Lfir and 
^HCO-i- should steepen from linear to super-linear over 
the range of galaxies presented in the CO panels. A care- 



ful, large-scale HC 0-|-(1 — ^ 0) survey similar to the work 
of Gao & Solomon (!2004ah on HCN(1 0) should reveal 
these trends. Lines with similarly low excitation temper- 
atures and intermediate critical densities like CS(1 — + 0) 
should behave analogously. 

4. DISCUSSION 

4.1. Implications for Kennicutt- Schmidt Laws and Star 
Formation Efficiencies 

Our results suggest that KS laws in different tracers 
naturally fall into two regimes, although there is a broad 
range of molecular tracers that are intermediate between 
the two extremes. Tracers for which the critical density is 
small compared to the median density in a galaxy repre- 
sent one limit. In these tracers, the light faithfully follows 
the mass, so the KS law measures a relationship between 
total mass and star formation. In any model in which 
star formation occurs at a roughly constant rate per dy- 
namical time, this must produce a KS law in which the 
star formation rate rises with density to a power of near 
1.5, and the ratio of star formation to luminosity rises 
as density to the 0.5 power. In terms of surface rather 

than volume densities, this implies S* oc S^^^/i^^/^. If 
we further add the observation that the scale heights h of 
the star-forming molecular layers of galaxies are roughly 
constant across ga laxy types, one form of the observed 
lKennic uT3 ('T998a!,fF) star formation law follows immedi- 
ately (Elmegreen 2002). More over, in a galactic disk , 
h (X '^g/n and n oc r^VQ (e.g. iThompson et al.l [20051) ; 
since in star- forming disks the Toomre-Q is about unity 
([Martin fc Kennicuttll200l[ ). substituting for h immedi- 
ately gives (X Efi, the alternate form of the [K"ennicuttl 
(1998 aibl) law. 

The other limit is tracers for which the critical density 
is large compared to the median galactic density. These 
tracers pick out a particular density independent of the 
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Fig. 4. — The same as Figure[3l but with constant V^^^l = 10* pc'^, and for the model parameters corresponding to "starburst" {dashed), 
"intermediate" (solid), and "normal" (dot-dashed) (see Table[T]l. Therefore, the middle solid line in each panel of Figure[3]is the same as 
the solid line in each panel in this Figure. 



mean or median density in the galaxy, and thus aU the 
regions they identify have the same dynamical time re- 
gardless of galactic environment. In this case the star 
formation rate will simply be proportional to the total 
mass of the observed regions, yielding a constant ratio of 
star formation rate to mass, as is obse rved for HON in 
the lo ca l universe (Fig Bl right panel; iGao fc SolomonI 

I2nn4allbl : Im et A\\\mm . 

We predict that there should be a transition between 
linear and super-linear scaling of L-pin with line luminos- 
ity at the point where galaxies transition from median 
densities that are smaller than the line critical density 
to median densities larger than the critical density. The 
HC0+(1 0) line, and other lines with similar critical 
densities, e.g. CS(1 0) and S0(1 — > 0), should show 
this behavior for galaxies in the local universe. The ob- 
served correlation between iHCO+ ^^'^ ^fir appears to 
be consistent with our prediction, although at present the 
data are not of sufficient quality to distinguish between 
a break and a single powerlaw relation. There are hints 
that the very highest luminosity star- forming galaxies, 
which all reside at high redshift and may well reach ISM 
densities not found in any local systems, show such a 
break in the IR-HCN correlation. 

One important point to emphasize in this analysis is 
that we have been able to explain the observed correla- 
tions between line and infrared luminosities, and hence 
between gas masses at various densities and star for- 
mation rates, without resorting to the hypothesis that 
the star formation process is fundamentally different in 
galaxies of different properties. Although uncertainties 
in both our model and the observations do not preclude 
an order-unity change in the star formation efficiency or 
SFRff as a function of Lfir, there is currently no evi- 
dence for such a change in the dat a , cont rary to claims 
made by, e.g. iGracia-Carpio et all ()2006f ). In fact, all 
of the observational trends are predicted by our sim- 
ple model with constant star formation efficiency. This 
is consistent with other lines of evidence that the frac- 
tion of mass at a given density that turns into stars 
is roughly 1% per free- fall time independent of density 
(iKrumholz &: Tanl[2007h . 

4.2. Does Star Formation Have a Fundamental Size or 
Density Scale? 



Based on the linear correlation between IICN(1 — > 0) 
luminosity and star formation rate, seen both in ex- 
ternal galaxies and in individual mo lecular clumps in 
the Milky Wav. lGao fc SolomonI (|2004a >b) and lWu et all 
propose that HCN(1 0) emission traces a fun- 
damental unit of star formation. They explain the linear 
IR-HCN correlation as a product of this; in their model, 
HCN luminosity correlates linearly with star formation 
rate because HCN luminosity simply counts the number 
of such units. 

Based on our analysis, we argue that this hypothesis is 
only partially correct. We concur with Gao & Solomon 
and Wu et al. that the HCN(1 — > 0) luminosity of a 
galaxy does simply reflect the mass of gas that is dense 
enough to excite the HCN(1 0) line. However, our 
analysis shows that this does not necessarily imply that 
this density represents a special density in the star for- 
mation process, or that objects traced by HCN(1 0) 
represent a physically distinct class. We show that a 
linear correlation between star formation rate and line 
luminosity is expected for any line with a critical density 
comparable to or larger than the median molecular cloud 
density in the galaxies used to define the correlation. It 
is possible that HCN(1 0)-emitting regions represent 
a physically distinct scale of star formation as IWu et all 
propose, but one can explain the linear IR-HCN correla- 
tion equally well if they are just part of the same contin- 
uous medium as the regions traced by C0(1 — > 0) and 
by other transitions. Even the star- forming clouds them- 
selves may simply be parts of a continuous distribution 
of ISM structures occupyin g the entire galaxy, as argued 
bv IWada fc Normlml ()2007[ ). In this case there need be 
no special density scales other than the mean and me- 
dian densities for the star-forming clouds on their largest 
scales, and the density at which star formation becomes 
rapid, converting the mass into stars in of order a free- 
fall time. This transition scale is unknown, but must 
be considerably larger than the density traced by HCN 
(jKrumholz fc T an 2007). 

4.3. Limitations and Cautions 
4.3.1. Self-Consistency 

As mentioned in § 13.21 our approach of leaving the gas 
temperature T and Mach number M as free parame- 
ters is not entirely consistent with our calculation of the 
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IR luminosity, since the IR luminosity and temperature 
are of course related. In principle, with a model of how 
the energy output from stars heats the dust and gas, to- 
gether with a structural model connecting the energy and 
momentum output from stars to the generation of turbu- 
lence, it should be possible to self-consistently compute 
both the gas temperature and the Mach number from the 
volum etric star formation rate (see, e.g. jThompson et al.l 
|2005[) . Such a model would return T and as a function 
of n and possibly other galaxy properties, while simulta- 
neously predicting a set of Kennicutt-Schmidt laws. 

If the line luminosity depended strongly on T or A^, 
or if one required knowledge of the temperature to com- 
pute the infrared luminosity of a galaxy, we would would 
have no alternative to constructing such a model if we 
wished to explain the observed IR-line luminosity cor- 
relation. However, we can avoid this by relying on the 
observationally-calibrated star formation-IR correlation, 
and because, as we show in Figures O [31 andlH the line 
luminosity varies quite weakly over a reasonable range of 
T and A4 for our chosen lines. For this reason, any model 
for computing T and as a function of galaxy proper- 
ties, if it were consistent with observations, would not 
significantly alter the IR-line luminosity correlation we 
derive. This is true, however, only for lines that require 
low temperatures to excite. As we discuss in § I4.3.2( 
lines that require higher temperatures to excite do de- 
pend sensitively on the temparature in the galaxy, and a 
model capable of predicting the IR-line luminosity cor- 
relation for these lines must also include a calculation of 
the temperature structure of the galaxy. 

4.3.2. Isothermality 

Our assumption of isothermality means that our anal- 
ysis will only apply to molecular lines for which the tem- 
perature Tup corresponding to the upper state energy is 
< 10 K, low enough to be excited even in the coolest 
molecular clouds in normal spiral galaxies. The reason 
for this is that at temperatures larger than T^p, the lu- 
minosity in a line generally varies at most linearly with 
the temperature. As the similarity between the results 
with our different galaxy models illustrates, changing the 
temperature within the range of ~ 10 — 50 K produces 
only a factor of a few change in the luminosity of the 
lines we have studied. In contrast, line luminosity re- 
sponds exponentially to temperature changes when the 
temperature is below the value corresponding to the up- 
per state's energy. This means that lines sensitive to 
high temperatures pick out primarily the regions that 
are warm enough for the line to be excited. Density has 
only a secondary effect. The emission will therefore re- 
flect the temperature distribution in star-forming clouds 
more than the density distribution, an effect that our 
isothermal assumption precludes us from treating. KS 
laws in high temperature tracers are likely to find lin- 
ear relationships between star formation rate and mass 
regardless of the critical density of the molecule in ques- 
tion because they will simply be correlating the mass 

of dust warmed to ^ 100 Kelvin, which is essentially 
what is measured by Tfir, with the mass of gas warmed 
to temperatures above Tup. However, our model will 
not apply to these lines, and for this reason we do not 
attempt to compare to observations using higher tran- 
sitions of CO (3 2, 4 ^ 3, 5 ^ 4, 6 ^ 5, and 



7-^6, whi ch have Tup = 33,55,83,116, and 154 K 
respe ctively; iGreve et al.ll2005L ISolomon fc Vanden Bout! 
[2005I) . CS(5 ^ 4) (Tup = 35 K: iPlume et al.lll997D . or 
other high temperature tracers. 

4.3.3. Molecular Abundances 

We have not considered density-dependent variations 
in molecule abundances. One potential source of vari- 
ation in molecular abundance is freeze-out onto grain 
surfaces at high d ensities and low temperatures (e.g. 
iTafalla et al.ll2004al lbl). Chemodynamical models suggest 
that freeze-out is not likely to become significant for ei- 
ther carbonaceous or nitrogenous species until densities 
10® cm~^ (jFlower et al.li2006. ), but may become se- 
vere at higher densities, so whether depletion is signif- 
icant depends on what fraction of the total luminosity 
would be contributed by gas of this density or higher 
were there no freeze-out. Figure [T] suggests that freeze- 
out is likely to modify the total galactic luminosity of 
CO, HCO+, and HON fairly little even at a mean ISM 
density of n = 10^ cm~'^, but may have significant effects 
for galaxies of larger mean densities or for lines for which 
the critical densities is comparable to the freeze-out den- 
sity. If freeze-out is significant, our conclusions will be 
modified. 

4.3.4. Atomic Gas 

In the simple model developed here, we have neglected 
the role of atomic gas entirely. Whether the density or 
surface density of atomic gas plays a role in controlling 
the star formation rate is subject to debate on both ob- 
servational and theor etical ground s feenn icutt 1998a, ij 
Wong fc Blitl 120021: liTever et all l2004l: iKom ugi et al] 
d iKrumholz fc McKedl2005l : iKennicutt et al.,.200';il . 
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so it is unclear how much a limitation this omission re- 
ally is. We can say with confidence that in molecule-rich 
galaxies, which prov i de almo st all the dynamic range of 
both the iKennicuttI ()1998al lH) correlation and the cor- 
relations illustrated in Figures [3] and IH the atomic gas 
plays almost no role simply because there is so little of 
it. Thus, our predictions should be quite robust, except 
perhaps at the very low luminosity ends of Figures [3] and 

m 

4.3.5. AGN Contributions 

A final point is not so much a limitation of our work as 
a cautionary note about comparing our model with ob- 
servations. We have included in our model IR luminosity 
only from star formation, and molecular line luminosity 
only from molecules in cold star-forming clouds. How- 
ever, an AGN may make a significant contribution to a 
galaxy's luminosity in the far infrared by direct heating of 
dust grains, and in molecular lines via an X-ray dissocia- 
tion region. Indeed, several of the systems with the high- 
est IR luminosities in Figures [3] and |4] are contaminated 
by AGN. As noted in §3.3) this complicates an assessment 
of our prediction of an up-turn in the T^^j^ — TpiR and 
^HCO+ ^ ^FiR correlations at high luminosity. This de- 
viation from linearity at high gas density (at fixed Vmoi, 
high Tfir) is an essential prediction of our model, but 
testing it relies on a careful separation of the contribu- 
tion of the AGN to b oth the I R and line luminosit ies (e.g. 
iMalonev et"al]|1996f ). In fact. lCarilh et al.l(|2005( ) discuss 
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the possibility that the AGN's contribution to the IR lu- 
minosity in these systems causes them to be above the 
local linear L'^q^ — Lfir correlation. Such a contamina- 
tion would mimic the prediction of our model. However, 
iGao et all (|2nol argue that the sub-millimeter galaxies 
in their sample are not AGN dominated and that just 
one of three quasars in their sample ( APM 08279 - H5255 ) 
has a large AGN IR component. See iGao et all (j2007D 
for more discussion. For these reasons we contend that 
although our model is consistent with the existing data, 
the current evidence for a break in the -^hcn ^ -^fir cor- 
relation should be viewed with caution and more data 
in high density/luminosity systems is clearly required to 
understand the role of AGN contamination in shaping 
the correlation. 

5. CONCLUSIONS 

We provide a simple model for understanding how 
Kennicutt-Schmidt laws, which relate the star formation 
rate to the mass or surface density of gas as inferred from 
some particular line, depend on the line chosen to define 
the correlation. We show that for a turbulent medium 
the luminosity per unit volume in a given line, provided 
that line can be excited at temperatures lower than the 
mean temperature in a galaxy's molecular clouds, in- 
creases faster than linearly with the density for molecules 
with critical densities larger than the median gas density. 
The star formation rate also rises super-linearly with the 
gas density, and the combination of these two effects pro- 
duces a close to linear correlation between star formation 
rate and line luminosity. In contrast, the line luminosity 
rises only linearly with density for lines with low critical 
densities, producing a correlation between star formation 
rate and line luminosity that is super-linear. 

Based on this analysis, we construct a model for the 
correlation between a galaxy's infrared luminosity and its 
luminosity in a particular molecular line. Our model is 
extremely simple, in that it relies on an observationally- 
calibrated IR-star formation rate correlation, it treats 
molecular clouds as having homogenous density and ve- 
locity distributions, temperatures, and chemical compo- 
sitions, and it only very crudely accounts for variations 
in molecular cloud properties across galaxies. Despite 
these approximations, the model naturally explains why 
some observed correlations between infrared luminosity 
and line luminosity in galaxies are linear, and some are 
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esis that the star formation process is somehow different, 
either more or less efficient, in different types of galaxies 
or for media of different densities. Instead, our model is 
able to explain the observed correlations using a simple, 
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One strong prediction of our model is that there should 
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powerlaw behavior. However, there is some preliminary 
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redshift galaxies more luminous than any found in the 
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vations it is difficult to rule out the alternative explana- 
tion of the break as arising due to a progressively rising 
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and at high redshift may be used to test our predictions 
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lines. 
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